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Robust PI Passivity-based Control of Nonlinear Systems: Application to 
Port-Hamiltonian Systems and Temperature Regulation 

S. Aranovskiy^ R. Ortega^ and R. Cisneros^ 


Abstract —This paper deals with the problem of control of 
partially known nonlinear systems, which have an open-loop 
stable equilibrium, but we would like to add a PI controller to 
regulate its behavior around another operating point. Our main 
contribution is the identification of a class of systems for which 
a globally stable PI can be designed knowing only the systems 
input matrix and measuring only the actuated coordinates. The 
construction of the PI is done invoking passivity theory. The 
difficulties encountered in the design of adaptive PI controllers 
with the existing theoretical tools are also discussed. As an 
illustration of the theory, we consider port-Hamiltonian systems 
and a class of thermal processes. 

I. INTRODUCTION 

In many practical applications the plant to be controlled 
has an open-loop stable equilibrium, e.g., at the origin, and 
we would like to add a controller to regulate its behavior 
around another operating point. In the case of linear sys¬ 
tems the dynamics remains invariant under coordinate shifts, 
therefore this task can be easily accomplished using the 
incremental model of the plant. Unfortunately, this is not the 
case for nonlinear systems, for which there is no obvious 
advantage of working with the incremental model. 

Another common requirement in applications is the use of 
standard proportional-integral (PI) controllers, which over¬ 
whelmingly dominate the industrial market. Although com¬ 
missioning a PI to operate around a single operating point is 
relatively easy, the performance will be below par in wide 
operating regimes, which is the scenario in modern high- 
performance applications. To overcome this drawback the 
current practice is to re-tune the gains of the PI controllers 
based on a linear model of the plant evaluated at various 
operating points, a procedure known as gain-scheduling. 
There are several disadvantages of gain-scheduling including 
the need to switch (or interpolate) the controller gains and 
the non-trivial definition of the regions in the plants state 
space where the switching takes place—^both problems are 
exacerbated if the dynamics of the plant is highly nonlinear. 
Another common commissioning procedure is to use auto¬ 
tuners, that heavily rely on the availability of a “good” linear 
approximation of the plant dynamics. To avoid the need to 
rely on linearization it is necessary to develop a procedure 
to design robust PI controllers for nonlinear systems with 
uncertain parameters. 
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Motivated by the discussion above in this paper we iden¬ 
tify a class of (input affine) nonlinear systems for which it is 
possible to design a PI controller with the following features. 
FI Regulation of the closed-loop system around the desired 
(non-zero) operating point should be guaranteed. 

F2 The PI controller should be robust, in the sense that re¬ 
duced knowledge of the system parameters is required. 
F3 To simplify the controllers commissioning, a well de¬ 
fined admissible range of variation for the PI pro¬ 
portional and integral gains, preserving closed-loop 
stability, should be provided. 

We propose the construction of a PI controller with the 
features F1-F3 for plants with unknown dynamics verifying 
the following assumptions. 

A1 The open-loop system is unknown but has a stable 
equilibrium at the origin. 

A2 The desired equilibrium belongs to the assignable set 
and admits a convex Lyapunov function. 

A3 The Lyapunov function is the sum of two functions, 
depending on the un-actuated and actuated coordinates, 
respectively. The first function is unknown while the 
second one is separable and linearly parameterized in 
terms of some unknown parameters. 

A4 The input matrix is constant, known and has n — m zero 
rows, where n and m are the dimensions of the state 
and input vectors, respectively. 

As indicated in the article’s title we exploit the fundamental 
property of passivity to design the proportional-integral pas¬ 
sivity based controller, which will be referred in the sequel 
as PI-PBC. The first step in the design is to, relying on A1 
above, invoke the celebrated theorem of Hill and Moylan 
[1] to identify a suitable passive output for the system, 
with storage function the Lyapunov function of the open- 
loop system. Since our interest is the regulation of non¬ 
zero equilibria, we then use the results of [2] to create a 
new passive output for the incremental model with a storage 
function that has a minimum at the desired equilibrium. As 
shown in [2], feeding back the passive output through a 
PI controller ensures stability of the desired equilibrium/or 
all positive definite PI gains. It is important to underscore 
that, since the passivity property has been established for the 
incremental model, the equilibrium can also be stabilized 
setting the control input equal to the (constant) value that 
assigns the equilibrium, say u*, which is univocally defined. 
However, this open-loop control will, obviously, be non- 
robust. In the robustness context of the present paper neither 
the plant dynamics nor the Lyapunov function are known 


and, consequently, we cannot compute neither the passive 
output nor u*. It is at this point that we invoke A3 and A4 
above to prove that, under these assumptions, it is possible 
to define ranges for the proportional and integral gains 
that make the PI-PBC implementable and, consequently, 
guarantee stability of the equilibrium. Another important 
feature of the proposed PI-PBC is that it requires only partial 
measurement of the state, namely, only the m state variables 
associated to the non-zero rows of the input matrix, referred 
in the sequel as actuated coordinates. 

Several practical applications of PI-PBC have been re¬ 
ported in the literature. This include, RLC circuits [3], power 
converters [4], fuel cells [5], electric drives [6] and mechani¬ 
cal systems [7]. In [8] a procedure to add an integral action to 
a non-passive output for a class of port-Hamiltonian systems 
was hrst proposed, and later extended in [9], [10]. To the 
best of our knowledge, the present paper is the hrst attempt 
to design PI-PBCs with guaranteed stability properties for 
systems with partially known dynamics. 

A natural question that arises at this point is the incorpo¬ 
ration of adaptation in the design of the PI (or PID). In the 
power converter application of [4] a parameter that enters in 
the dehnition of the passive output, i.e., the load resistance, 
is adaptively identihed—however, all other parameters are 
assumed to be known. In the interesting paper [11] it is 
shown that it is possible to adaptively estimate u* for a 
general nonlinear system with scalar input, keeping the 
estimate in a known interval, provided the passive output is 
known. In spite of a large number of publications the problem 
of designing a provably stable adaptive PID for systems with 
unknown parameters remains, as far as we know, open. The 
difficulty of this task was identihed already in 1984 in [12]. 
As is well-known [13], the stability of indirect adaptive 
methods relies on parameter convergence that, in its turn, 
requires persistency of excitation—a property that is not 
satished in the regulation tasks where PI control is effective. 
On the other hand, the application of direct methods is 
stymied by the absence of a suitable parameterization of 
this structure-constrained controller. For the PI-PBC studied 
in this paper the main difficulty is the need to estimate 
two objects, that appear multiplicatively in the Lyapunov 
analysis: the passive output and the ideal control signal u*. 
This point is further elaborated in Subsection IV-BI 

The remaining of the paper is organized as follows. 
Section [B] presents the problem formulation. To streamline 
the presentation of the main result, which is given in Section 
IIVI some preliminary lemmata are given in Section [III] In 
Section lY] we discuss the reasons that stymie the use of 
adaptation and the inability to state a robustness result based 
on continuity and approximate prior knowledge of the plant. 
Section |Vl] is devoted to application of the proposed PI-PBC 
for port-Hamiltonian (pH) systems [1] and a temperature 
regulation problem. The paper is wrapped-up with conclud¬ 
ing remarks in Section IVIII 


zeros. Given £ K, i € n := {1, ..., n}, we denote with 
col(ai) the n-dimensional column vector and diagjai} the 
diagonal nxn matrix with elements a^. For x £ R", |a;| is the 
Euclidean norm. For mappings of scalar argument p : R 
R®, g' and g" denote first and second order differentiation, 
respectively. For mappings / : R” —>■ R, V/ := and 

V^/ := 1^. For the distinguished element x* £ R" and 
any mapping F : R” —>• R® we denote F* := F{x*) and the 
error signal F{x) := F{x) — F*. 

II. PROBLEM FORMULATION 

In this section we formulate the control problem addressed 
in the paper, enunciate the assumptions made on the plant to 
solve it and make some remarks on these assumptions. 


A. Robust PI control problem 

Consider the nonlinear, input affine, system 

X = f{x) + Gu, (1) 


where x £ R", u £ R™, n > m, / : R" —>• R" is an 
unknown smooth mapping, G £ is constant verifying 

rank(G) = m. 

The following is a key assumption made throughout the 
paper. 

Assumption 1: The matrix G has n — m zero rows. With¬ 
out lost of generalit}0 it is assumed of the form 

G= , (2) 

G2 


where G 2 £ R™xm icnown. 

This assumption can be easily obviated introducing state 
and input changes of coordinates. Indeed, it is well-known— 
see, e.g.. Theorem 2 of Section 2.7 of [14]—that for any full 
rank, matrix G £ K"xm (jjgj-g exists (elementary) full rank 
matrices T £ R"^" and S £ 

_ ll(n—m)xm 

Consequently, introducing z = Tx and v = S~^u the system 
O takes the desired form 


i = w{z) + 


(n—m) xm 


V, 


where w{z) = Tf{T~^z). We should note, however, that 
a change of state representation destroys—in general—the 
original structure of the system, whose knowledge may 
be critical for the verification of the second assumption 
below. This fact is clearly illustrated in the physical examples 
considered in Section |Vl] For this reason, we prefer to leave 
it as an standing assumption. 

Motivated by Assumption[T]we find convenient to define a 
partition of the state vector into its un-actuated and actuated 


Notations /„ is the nXn identity matrix and 0„xs is an nxs subsection and Subsection |Wa] for more general 

matrix of zeros, 0„ is an n-dimensional column vector of forms of G. 
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C. Discussion 


X = 

'^u 

1 ■ — 

Xi 

X2 

, Xa : = 

^n—m+1 

^n—m+2 




^n—m 


Xn 


It is assumed that only Xa is available for measurement. 

The unforced system, that is, x = f{x), has a stable 
equilibrium at the origin with a partially known Lyapunov 
function. We are interested in controlling the system with 
a PI at a non-zero equilibrium—a situation that arises in 
most practical applications. Thus, we are given a desired 
equilibrium point, x* G ffi.", and the control goal is to ensure 
stability of this equilibrium using a PI control law of the form 

i = -Ki'lp{Xa,X*) 

U = -Kp-4){Xa,X*) + Z, 

where z G K™ is the controller state, Kp G and 

Kj G are tuning gains and ■0 : K™ x R" —> R™ 

is a mapping designed with the partial knowledge of the 
aforementioned Lyapunov function. 

The following, practically reasonable, assumption is made 
throughout the paper. 

Assumption 2: The desired equilibrium point x* belongs 
to the assignable equilibrium set, that is, 

x* G f := {x G R" I [ In-m I ] /(x) = 0} . 

( 3 ) 

B. Assumptions on the open-loop plant 

The following assumption identifies the class of vector 
fields f{x) for which we provide an answer to the problem. 

Assumption 3: For the system O there exists a twice- 
differentiable, positive definite function H : R" -G R>o, 
verifying the following. 

(i) [VH{x)]^f(x) < 0. 

(ii) [VH{x) - yH{x*)]^fix) =: -Q{x) < 0. 

(iii) The function H{x) is of the form 

H{x) = Hu[Xu) + Ha{Xa) 


The following remarks regarding the assumptions are in 
order. 

R1 Although the set 8. is not known, it is reasonable to 
assume that we have enough prior knowledge about the 
plant to select the desired operating point as a feasible 
equilibrium. Hence, Assumption |2] is reasonable. 

R2 A corollary of Assumption |2] is that the constant input 
u*, that assigns the equilibrium, is uniquely defined as 

:= (GJG2)”'[O^x(n-m) Gl]f\ (5) 

Notice that, without knowledge of f{x), this constant 
cannot be computed. 

R3 Since the open-loop system ([TJ has a stable equilibrium 
at the origin Assumption [3 (i) follows as a corollary 
of Lyapunov’s converse theorems [15]. As will become 
clear below Assumption |3] (ii) and (iv) are required to 
prove passivity of the incremental model as done in [2]. 
R4 We underscore that no assumption, beyond twice dif¬ 
ferentiability and convexity, is imposed on the un¬ 
known component Hu{xu) of the Lyapunov function 
of the open-loop system H(x). On the other hand, 
stricter conditions are imposed on the second compo¬ 
nent Ha{xa), with uncertainty captured by the unknown 
constants di. 

R5 Assumptions [3] (iii) and Assumption [T] are the key 
requirements imposed on the plant to design the robust 
PI-PBC. This assumption is satisfied by a large class of 
physical systems, including a class of port-Hamiltonian 
[1] and thermal systems studied in Section |VT] 

R6 Regarding Assumptions [T] in the more general case 
when G is not of the form dU an additional shuffling of 
the rows of G is needed in the design. This procedure 
is explained in Subsection IV-AI 
R7 For quadratic Lyapunov functions of the form H(x) = 
x^Px, with P > 0, Assumption |3] (ii) is satisfied if 
the open-loop system is convergent in the sense of 
Demidovich [16]. That is, if it satisfies 

PV/(x) + [V/(x)]^P < 0. 

III. PRELIMINARY LEMMATA 

Unless otherwise indicated, throughout the rest of the 
paper Assumption [T] holds. Define for the system ([T]l the 
output 

y = G^yHix) = GjD^Xa), (6) 


with 

n 

i=n—m+l 


where 


( 4 ) 


D 


dn—m+l 

0 


0 

dn—m+2 


0 

0 


> 0 


where the function Hu : R"”™ —R and the constants 
di > 0 are unknown but the the functions (/)i : R —5> R 
are known. 

‘^{Xa) 

(iv) The functions Hu{xu) and (j)i{xi) are convex. 


0 0 ... dn 

4’n-m+l{Xn-m+l) 

4>'iXn) 














A corollary of the theorem of Hill and Moylan [1] is that, if 
Assumption [ 3 ] (i) holds, the system ([T]), (|6]) defines a passive 
mapping y with storage function H{x). 

To operate the system at a non-zero equilibrium it is 
necessary to shift the minimum of the storage function and 
define the passivity property between the incremental input 
and the output error. Towards this end, we recall Proposition 
1 of [2] and state it as a lemma below. 

Lemma 1: Consider the incremental model of the system 

O, ® 

X = f{x) + Gu* + Gu, 

-r ~ ( 7 ) 

e = GjD^Xa), 

where u = u — m* is the incremental input. Under Assump¬ 
tions [T]-[3]the mapping u 1 —>■ e is passive with storage function 
[/ : M” ^ ]R>o given by 

U{x) = H{x)-xlVH:-xjD^^ + k, ( 8 ) 


where A: is a constant that ensures U{0) = 0. More precisely, 

if = —Q(x) + e^u, (9) 

where Q(x) is defined in Assumption [3 (ii). 

One of the main interests of passive systems is that they 
can be globally stabilized with PI controls (with arbitrary 
positive definite gains). This well-known fact is stated in 
the lemma below, whose proof is given to streamline the 
presentation of our main result. 

Lemma 2: Consider the system O verifying Assumptions 
[TJ^in closed-loop with the PI-PBC 

e = G^D^^Xa) 

z = -Kie (10) 

u = —Kpe + z. 


For all positive definite gain matrices Kp G 
AT/ G all trajectories are bounded, the equilibrium 

point (x,z) = (x*,u*) is globally stable (in the sense of 
Lyapunov) and the augmented error signal 


• — 


Q{x) 

e 


( 11 ) 


where Q{x) is defined in Assumption [3 (ii), verifies 


lim ea{t) = 0. (12) 

t—^CO 


Moreover, if Ca is a detectable output for the closed-loop 
system then the equilibrium point is asymptotically stable. 

Proof: Defining z := z — u* the last two equations of 
the controller (fTOl l may be written in the form 


z = —Kje 
u = —Kpe + z. 


( 13 ) 


Consider the Lyapunov function candidate 

W{z,x) = U{x) + ^2^Aiz, (14) 


where A/ > 0. The time derivative of the Lyapunov function 
along the trajectories of the closed-loop system is 

W = —Q{x) + u + z^ Ajz 


= —Q{x) — Kpe + z^e — z^ KiKje. 


( 15 ) 


Setting A/ = ATj ^ yields 


W = —Q{x) — e^ Kpe. 


The proof is complete invoking standard Lyapunov argu¬ 
ments [15]. ■ 


IV. THE ROBUST PI-PBC 

As indicated in R4 of Subsection III-CI the matrix D is 
unknown. Hence, the error signal e cannot be constructed 
and the PI-PBC (ITOl) is not implementable. This motivates 
our main result given below. 

Proposition 1: Consider system O verifying Assump¬ 
tions 1-3 in closed-loop with the robust PI-PBC 

U = -Kp^lXa) + Z 

~ (16) 

Z = -Ki^{Xa), 

with the controller gains 

Kp = G2 p 

Ki = Gf^Ti. (17) 

For all diagonal, positive definite matrices Lp S and 

L/ e ]^mxm Jj^yg Jjjg following. 

(i) All trajectories are bounded and the equilibrium point 
{x,z) = {x*,u*) is globally stable (in the sense of 
Lyapunov). 

(ii) The augmented error signal Cq defined in (ITTI) verifies 

(Hill. 

(iii) If Ca is a detectable output for the closed-loop system 
then the equilibrium point is globally asymptotically 
stable. 

Proof: Some simple manipulations prove that 

Kp^xa) = Gf^TpD-^Gf^Gj D^xa) = Ape, (18) 

where we defined the matrix 

Ap :=Gf^rpD-^Gf~^, (19) 

and used the definition of e in ([TOb . Invoking Sylvester’s Law 
of Inertia [14], and the fact that Lp and D are diagonal and 
positive definite, we have that Ap > 0. 

Next choose 

Ai ■.= GjDrj^G2, ( 20 ) 

that is, also, positive definite for all diagonal, positive definite 
L/. Then 

AiKi^ix) = GjD^{xa) = e. (21) 

Replacing dTSl ) and (ISTT i in the controller equations yields 

u = —Ape -\- z 
z = —Af^e. 






Consequently, the time derivative of the Lyapunov function 
O becomes now 

W =-Q{x) - kpe, (22) 

completing the proof. ■ 

To obtain an implementable version of the robust PI- 
PBC it was necessary to carry-out two tasks. First, to make 
the damping injection, introduced by the proportional term, 
function of the unknown matrix D. Indeed, replacing (fT9l l in 
(flSl l we get 

Kp^{x) = G'^^TpD-^G^^e. 

Second, make the gain A/ of the Lyapunov function (fl4l i 
also a function of D —see (EOll. 

An important observation is that, even though the con¬ 
troller only requires measurement of the actuated terms of 
the state Xa, it achieves regulation of the full state vector. 


V. ADDITIONAL REMARKS ON THE PI-PBC 


In this section we explain how to proceed when G is not 
of the form ©, discuss the reasons that stymie the use of 
adaptation and the inability to state a robustness result based 
on continuity and approximate prior knowledge of the matrix 
D. 

A. General G (with n — m zero rows) 


Instrumental to design the robust PI-PBC was the partic¬ 
ular form of H{x) defined in Assumption [3 (iii). In view 
of the construction of the robust PI-PBC, it is clear that if 
G is not of the form ^ the assumption must be modified 
redefining the actuated and un-actuated coordinates. 

To avoid cluttering the notation we will explain the proce¬ 
dure only for the case when n = 3 and m = 2—the general 
case follows verbatim. Assume, furthermore, that G is of the 
form 


G = 


9 I 


93 


The form of H{x) given in Assumption [3 (iii) must be, 
accordingly, modihed to 


H{x) = Hu{X 2 ) + di(j){xx) + d 3 (j){X 3 ). 


In this case the passive output e for the incremental model 
becomes 

G^[VH{x)-VH(x*)] = Gs ° 

where 

Gs-=[ 91 I 93 ]■ 

The robust PI-PBC is given by 

§i(a;i) 

$3(2^3) 

$i(xi) 

‘k3(2;3) 


U = -G-^Tp 

z=-Gj^Ti 


$i(xi) 

^’3{2;3) 


where Ep and E/ are arbitrary, diagonal, positive definite 
matrices. 

Before closing this subsection we remark that our con¬ 
struction critically relies on the assumption of existence of 
n — m zero rows in G. Indeed, it is possible to show that if 
this is not the case, even assuming H{x) of the form 

n 

H{x) = y^^di(j)i{xi) 

i=l 

and defining = diagjdi}, it is not possible to find an 
mxm positive definite matrix A, which will depend on Z?„, 
such that the matrix KG^Dn is independent of £>„. The 
fact that this is not possible for all matrices G is obvious 
considering the counterexample G = col(l,l). Hence, the 
assumption of existence of n—m zero rows in G is necessary 
to solve the problem. 

B. Difficulties for adaptation 

A natural alternative to the robust PI-PBC presented above 
is to assume a parametrisation of f{x) and try to estimate 
this parameters or, in a direct approach, estimate the matrix 
D that defines the passive output. The indirect approach, as 
is well-known, relies on parameter convergence that requires 
persistency of excitation—a property that is not satisfied in 
the regulation tasks where PI control is effective. 

Let us see what are the difficulties for the application of 
a direct adaptation approach. Towards this end, we propose 
the adaptive PI-PBC 

D = F{x,z) 
e = GjD^Xa) 
z = —Kje 
u = —Kpe + z, 

where the parameter adaptation law F : R" x R™ ^ j^mxm 
is to be definedll Defining e := e — e the last two equations 
of the controller may be written in the form 

z = —Ki{e + e) 
u = —Kp(e + e) + z. 

The time derivative of the Lyapunov function (fT4li with A/ = 
is now 

W = —Q{x) — Kpe — u^e 

= —Q(x) — Kpe — G^ D^{xa) 

where we underscore the presence of the last right hand term. 
If u were known the standard procedure of augmenting the 
Lyapunov function with a term trace(II^L)) and cancelling 
the sign-indefinite term with a suitable choice of F{x,z) 
would do the job. Alas, u* is not known, hence this approach 
is not feasible. 

Adding an adaptation for the constant u* is also not 
a trivial task, because of the bilinear nature of the joint 
estimation problem. 

^Notice that, in contrast to the robust PI-PBC, we have assumed that the 
full state is measurable. 












C. Comments on robustness based on continuity 

The availability of a bona fide Lyapunov function for 
the known parameters Pl-PBC, i.e., W{x,z), suggests that 
stability will be preserved if the matrix D is replaced by a 
“good”, constant estimate of it, say Dq. More precisely, it is 
expected that replacing the controller (fTOl i by 

eo = GjDo^{xa) 

Z = -KiCq 
u = -KpCo + z, 

where 

D = Do + A, A := diag{(5J 

would ensure stability if |col((5i)| is sufficiently small. Un¬ 
fortunately, since the Lyapunov function is not strict, this 
conjecture cannot be analytically validated. Indeed, in this 
case the time derivative of the Lyapunov function (fT4li with 
A/ = is now 

W = —Q{x) + u — z^(e — GJ A^(xa)) 

= -Q(x) - e^Kpeo - {Kpco - z)^GJA^{xa). 

While the term Cq KpGjA^{xa) can be dominated for 
“small” A, there is no way we can dominate the remaining 
term z^KjGj A$(a;a) and the Lyapunov analysis cannot be 
completed with standard techniques. 

This unfortunate situation does not mean, of course, that a 
continuity result of this type cannot be established. It simply 
reveals our inability to do it with the tools used to analyze 
the ideal case. 

VI. EXAMPLES 

Two examples of physical systems which are amenable for 
robust PI-PBC are given below. Attention is concentrated on 
the verification of Assumption [3 Hence, unless otherwise 
indicated. Assumption [T] is not imposed. 

A. A class of port-Hamiltonian system 
Proposition 2: The pH system 

X = {J - n)\/H{x) + Gu (23) 

with constant interconnection J = —J^ and damping 72. = 
127 > 0 matrices satisfies Assumption [3 (i) and (ii). 

Proof: Assumption [3 (i) is obviously satisfied because 

V^H{x){J - n)VH{x) = H{x)n\/H{x) < 0. 

Similarly, Assumption [3 (ii) holds since 

[SJH{x) - VH{x*)\^{J - TlfiVHix) - VH{x*)] = 
-[VH{x) - VH{x*)]^n[7H{x) - VH{x*)] < 0. 

■ 

B. Application to temperature regulation 

In this subsection we design a robust PI-PBC for the 
temperature regulation of a class of thermal systems—the 
so-called, rapid thermal processes. 


1) System Description: Similarly to [17], [18] we con¬ 
sider the following model of rapid thermal processes 

T = Ai [^(T) - ^{Trad)] + A 2 (T - + Gu, (24) 

where T G K.>g represents the vector of temperatures, 
T'(T) := col(r/) and Trad, 7onv G R>o are, respectively, 
the vectors of temperatures related to the radiation heat 
emission from environment and the convection air flows. The 
constant matrices Ai, A 2 G are Hurwitz and contain 

the radiation and the convection heat transfer coefficients. 
Also, the entries of G € correspond to the heat 

transfer coefficients of the heating elements. Einally, u G R"* 
is the controlled power applied to the heating elements. 

In the model above, as in [18], it is considered that (l24l) 
is heated almost uniformly so that the contribution of energy 
from conduction is too small with respect to the radiation 
transfer. Hence, the conduction heat transfer is neglected. 

To simplify the notation we re-write the system (i24l i in 
the form 

f = AiT'(T) -f A 2 T + E + Gu, 

where 

E := -AiA>{Trad) - A2Taouv 

Unlike Ai, A 2 and E that are highly uncertain, the input 
matrix G —that is defined by the induced heat flow—can be 
precisely identified. The control objective is then to design 
a robust PI, i.e., that does not require the knowledge of the 
uncertain parameters, to regulate the process around some 
desired temperature, which is not equal to the open-loop 
equilibrium, but belongs to the assignable equilibrium set, 
that is, 

T* g{T G R|o I G-^[AiT'(T) -f A 2 T -p l;] = 0} , (25) 

where G-^ G ^{n-m)xn ^ full-rank left-annihilator of G. 

To place the problem in the context of Proposition 1 we 
first shift the equilibrium of the open-loop system to the 
origin and then proceed to verify Assumption [3 Eor, we 
introduce the standard change of coordinates 

x = T-T, 

where T is the open-loop equilibrium that satisfies 

AiAi{T)+A2f + E = 0. (26) 

Thus, the system (i24l i in the new coordinates takes the form 
([T]i with 

f{x):=AiA>{x + T) + A 2 {x + f) + E, (27) 

Associated to the desired temperature T* we define the 
equilibrium to be stabilised 

X* := T* - T. 


( 28 ) 


2) Passivity of the thermal system: The lemma below 
identifies conditions under which the system (l24l i satisfies 
Assumption [3 without imposing Assumption [T1 that is, 
avoiding the partition of the coordinates into actuated and 
un-actuated. Towards this end, the following assumption is 
needed. 

Assumption 4: The matrix Ai is diagonally stable [19]. 
That is, there exists P € P = diagjpi} > 0 such that 

PAi+AjP=:-2S<0. (29) 

Moreover, the matrix A 2 verifies 

AjPdiag{Tf} + dmg{Tf}PA 2 < 0. (30) 

Conditions for diagonal stability of a matrix have been 
studied intensively, see [19] for a survey. Necessary and 
sufficient conditions were first reported in [20]—see also [21] 
for a simple proof. A sufficient condition, given in [22], is 
that it is a Metzler matrix (namely, its non diagonal elements 
are nonnegative). 

Since A 2 is Hurwitz, condition dSOl l is trivially satisfied if 
A 2 is diagonal, which is the case in some physical examples. 

Lemma 3: If Assumption |4] holds the vector field (IZTl i 
satisfies Assumption |3] with 

n 

H{x)=Y, Pi4>i{xi) + k, (31) 

where ^ 

+ T^)^ - ^fi{fi)x^ (32) 

5 

and 

i=l 

Proof: Point (iii) of Assumption [3] is trivially satisfied 
by dlD. 

We proceed now to prove point (i). Replacing ( l32b in ( |3T] ) 
and grouping terms yields 

1 ” 

H{x) = Pi{xi + Ti)^ - x^P'1!{T) + k, 

Now, notice that 

VH{x) = P^{x), 

where 

$(a:) := T'(a; + T') - T'(T). (33) 

On the other hand, from (i26T l it follows that the systems 
vector field may be written as 

f{x) = AiT>(a;) + A 2 X. 

Consequently, 

\VH{xyy f{x,9) = {x)P\Ai^{x) + A 2 x\ 

= —+ $^(a;)PA2a:, 

where we have used (i29] l to obtain the second identity. Now, 
condition (O ensures that the function : R" —>• R" 

h{x) := AjPT'(a;), 


is monotonically decreasing [16]. That is, for all a, 6 € R", 
[h{a) — h{b)y {a — b) < 0. 

Consequently, 

^^{x)PA 2 X = [hix + f)- h(T)yx < 0 

completing the proof of point (i). 

To prove point (ii) we notice that 

f{x) = Ai$(a;) + A 2 X, 

while 

VH{x) - VH{x*) = P^{x). 

Hence, the claim is established invoking the same arguments 
used above and defining 

Q(x) = ^^{x)S^{x). 

Finally, the second derivative of the functions 4>i{xi) yields 
f'l{x,) = d{x,+%f =dTf, 

which is non-negative because Ti > 0. Hence, the functions 
fiixi) are convex as requested by condition (iv) of Assump¬ 
tion 13 This completes the proof. 

■ 

Direct application of Lemma 1 leads to the following. 
Corollary 1: If Assumption |4] holds, the thermal system 
dull defines a passive map u e with storage function 
U{x), where 

e = G^P$(x) 

U{x) = H{x)-x^P^{x*)-H{x*) + {x*y P^x*). 

3) Robust PI-PBC of the thermal system: To present 
the robust PI-PBC for systems verifying Assumption [T] we 
partition the vector of temperatures into its un-actuated and 
actuated components 





Ti 


2n—m+1 

T = 

■ ■ 

. Ta . 

T — 

1 U • — 

T 2 

, Ta-.= 

m+2 




Tn-m 




partition P as 

p _ ®(n—m)xm 

^ n r) ’ 

y^mx{n—m) ^ 

and do the same with the vector function T'(T). 

The following proposition is a consequence of Lemma |3 
and Proposition [T] 

Proposition 3: Consider the system dull verifying As¬ 
sumptions [Hand |4l Fix any desired temperature T* verifying 
d25l l and define the PI-PBC 

U = -Kp^aiTa) + Z 

Z=-Ki^aiTa), 

and the controller gains Kp and Kj are given by dnll- For 
all diagonal, positive definite matrices Fp G 
F/ G ]^mxm trajectories are bounded and the equilibrium 
point {T,z) = {T*,u*) is globally asymptotically stable. 










Proof: The proof of stability is established invoking 
item (i) of Proposition 1 and identifying 

^a{Xa)\xa^=Ta-Ta ~ ^a{Ta)- 

To prove asymptotic stability we invoke item (ii) and observe 
that the augmented error signal (fTTTi is given in this case by 

^ \^^{T)S 
“ [ GjD 

Since Cq verifies (fT^ and S is positive definite we conclude 
that 'T(T(f)) —0 and consequently T{t) —>■ T*. ■ 

Physically, considering matrix G as (l2]l means that for m 
heating elements there are n — m measured points that are 
not directly heated by these elements. 

VII. Concluding Remarks 
In this work we identify a class of nonlinear systems for 
which it is possible to design robust PI controllers with 
guaranteed stability properties. The class consists of input 
affine systems with known, constant input matrix G and 
n — m zero rows. We assume that only the states associated 
to the non-zero rows of G are measurable.The systems have 
an open-loop stable equilibrium, but is different from the 
desired operating point. To handle this situation, we follow 
[2] and generate new passive outputs for the incremental 
model, hence the name PI-PBC. Associated to the open- 
loop stable equilibrium a Lyapunov function of the form (01 
is assumed to exist. We underscore that, besides convexity, 
there is no assumption on the function Hu{xu), which is 
unknown. Moreover, the controller does not require the mea¬ 
surement of Xu- The functions (j)i{xi) are assumed convex 
and known, but the coefficient di are unknown. Under these 
conditions, we show that, for a well identified class of PI 
tuning gains, see dnii, global stability of the proposed PI- 
PBC is guaranteed. Conditions that ensure global asymptotic 
stability, are also derived. 
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